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Abstract 

The impact of recent precision measurements of DIS structure functions and inclu- 
sive jet production at the Tevatron on the global QCD analysis of parton distribution 
functions is studied in detail. Particular emphasis is placed on exploring the range 
of variation of the gluon distribution G(x,Q) allowed by these new data. The strong 
coupling of G(x, Q) with a s is fully taken into account. A new generation of CTEQ 
parton distributions, CTEQ4, is presented. It consists of the three standard sets (MS, 
DIS and leading order), a series that gives a range of parton distributions with cor- 
responding a s 's, and a set with a low starting value of Q. Previously obtained gluon 
distributions that are consistent with the high E t jet cross-section are also discussed 
in the context of this new global analysis. 
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1 Introduction 



Lepton-lepton, lepton-hadron, and hadron-hadron interactions probe different and comple- 
mentary aspects of Quantum Chromodynamics. Each of these interactions provides a win- 
dow onto the elementary interactions of quarks and gluons and probes the running coupling 
and quark masses. The corresponding calculations are performed using perturbative Quan- 
tum Chromodynamics (pQCD). In addition these interactions probe the partonic structure 
of hadrons, as represented by the parton distribution and fragmentation functions |l] ||. 
These functions are essentially non-perturbative. There are, of course, large areas of overlap 
between various processes, which provide impressive consistency checks of the theory. In 
the first approximation, lepton-lepton processes provide clean measurements of basic pa- 
rameters such as quark charges, the strong coupling a s (Q), and fragmentation functions 
of partons into hadrons. Deep inelastic scattering structure functions and lepton-pair pro- 
duction cross-sections in hadron collisions provide the main source of information on the 
quark distributions f q (x,Q) inside hadrons. At leading order, the gluon distribution func- 
tion G(x, Q) enters directly in hadron-hadron scattering processes with direct photon and jet 
final states. In a global QCD analysis incorporating all these processes, one tries to exploit 
the strengths of each process in a uniform framework. Modern analyses are carried out to at 
least next-to-leading order (NLO), thus {a s (Q), f q (x,Q), G(x,Q)} all contribute and mix 
in the theoretical formulas for each process. However, the broad picture outlined above does 
reflect the main roles the various processes play in the analysis. 

Direct photon production has long been regarded as potentially the most useful source 
of information on G(x, Q). Fixed-target direct photon data, especially those from WA70 B, 
have been widely used in existing global analyses. However, there are a number of theoreti- 
cal uncertainties which affect the predictions of the normalization and slope of the measured 
direct photon p t spectrum. These effects include: (1) the sensitivity of the theoretical calcu- 
lations to the choice of factorization and renormalization scales [|J]; (2) k t broadening of the 
initial state partons due to soft gluon radiation [|J; and (3) photon fragmentation uncertain- 
ties || and the related issue of photon isolation cuts 0. When all these uncertainties are 
taken into account, existing direct photon data do not place as tight constraints on the gluon 
distribution as is commonly believed 0. Full exploitation of the potential of this process 
in a QCD global analysis will require significant progress in the understanding of the above 
issues. 

An important process that is sensitive to gluons is jet production in hadron-hadron col- 
lisions. In leading order, the cross-section is proportional to a 2 s {Q) G(x, Q) G(x', Q) and 
a 2 s {Q) G(x, Q) q(x', Q) for the gluon-gluon and gluon-quark scattering subprocesses respec- 
tively. Experimental measurement of various inclusive jet cross-sections has progressed to 
an increasingly quantitative level in recent years. For instance, at the Tevatron, good data 
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on single jet production are now available over a wide range of transverse energy, 15 GeV 
< E t < 450 GeV |7|, ||. NLO QCD calculations of jet cross-sections have also reached 
a mature stage || [10], |H| . Many issues relating to jet definition (which is important for 
comparing theory with experiment) encountered in earlier stages of jet analysis have been 
extensively studied and are better understood. For the moderate to large E t range, the scale 
dependence of the NLO inclusive jet cross section turns out to be relatively small | 12|| . Thus, 
it is natural that inclusive jet data should now be incorporated into a global QCD analysis, 
and that these data should play a role in constraining the gluon distribution G(x, Q). 

We have carried out a first systematic study of this problem using the CTEQ global anal- 
ysis framework |T3|]{] In this paper, we shall focus on the question: How well can the gluon 
distribution be determined as the result of recent advances in experimental measurements? 
We discuss phenomenological issues pertinent to extracting G(x, Q) in the global analyses. 
These factors are systematically taken into account in a series of analyses to gain insight on 
the current range of uncertainties on G(x, Q). We found that recent, more precise, DIS data 
have a significant influence in narrowing down the parton distribution functions (PDF's), 
including G(x, Q); and the inclusion of inclusive jet data from hadron colliders further solid- 
ifies knowledge on G(x, Q) over a wide range of x. As the result of this study, we present 
new sets of CTEQ parton distributions (in MS, DIS and LO schemes) as well as a series of 
distributions which give a range of variation of PDF's consistent with current data. We give 
quantitative information on how these distributions compare to the data sets used in the 
analysis. In addition, we provide a new set of PDF's with a low initial Qq = 0.5 GeV 2 ; and 
we discuss the previously obtained gluon distributions designed to accommodate the high 
E t jets H in the context of the CTEQ4 analysis. 



2 Issues on the determination of the gluon distribution 

In pQCD, the gluon distribution function is always accompanied by a factor of the strong 
coupling (i.e. it appears as a s G(x, Q)), both in the hard cross-sections and in the evolution 
equation for the parton distributions. Thus, the determination of a s and G(x, Q) is in 
general a strongly coupled problem. In principle, a s can be independently extracted from 
e + e~ collisions, or in sum rule measurements in deep inelastic scattering. G(x, Q) can then 
be determined in a global analysis, along with the quark distributions f l (x,Q), by treating 
a s as known. Alternatively, one can try to determine a s , G(x, Q) and the quark distributions 
at once in a global analysis. This relies on the full (x, Q) dependence of the wide range of 
data to differentiate a s (which controls the overall Q-dependence of all quantities) from the 
parton distributions (which depend on both x and Q). This method is not as "clean" as the 

X A similar analysis was carried out earlier 0, focusing on the interpretation of the "high pt excess" seen 
in the CDF jet measurement @. See Sec. ^ for a discussion of the relation of [|| to the current analysis. 
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first approach, and it will not become precise until the global analysis system has become 
better constrained. Eventually, however, it is important to demonstrate that the same value 
of a s consistently describes all the processes included in the global analysis. Hence, the two 
approaches are indeed complementary. 

It is well known that, at present, the value of a s determined at high energy colliders, 
especially LEP, is generally higher than that obtained from analyses of fixed-target DIS data 
IIIH . Since global QCD analyses are up to now dominated by the copious high statistics 
DIS data, they favor values of a s close to the lower "DIS value" . This situation may change 
when more and more quantitative results from hadron collider processes, such as inclusive jet 
and direct photon production, are included in the global analysis. In the following, we shall 
explore the range of variation of G(x, Q) when the value of a s is varied within the currently 



accepted region, which we shall take to be 0.105 < a s (Mz) < 0.122. [15] The problem of 
the determination of a s in global analysis and the question about consistency of a s among 
different processes will be considered in a subsequent study [[U| . 



For a quantitative study of G(x,Q), another relevant consideration is: How does the 
choice of parametrization of the initial gluon distribution at some Q = Qo affect the results? 
All global analyses use a generic form: 

G(x, Qo) = A x M (1 - x) M P(x; A 3 , ..) (1) 

with Ax 2 being physically associated with small-x Regge behavior and large-x valence count- 
ing rules respectively; and P(x; A 3 , ...) being a suitably chosen smooth function depending 
on one or more parameters. In general, both the number of free parameters and the func- 
tional form can have an influence on the global fit. In the CTEQ3 analysis an effort was 
made to minimize the number of free parameters, resulting in an economical set whereby 
Af = Al ea , and Pcteqz{ x 'i A 3 ) = 1 + A 3 x. We shall refer to this choice as the minimal set 
in the following discussions. In the literature, more degrees of freedom have been assigned 
to G(x, Qo). For instance, in CTEQ2 ]13[ and in recent MRS fits Af is allowed to vary 



independently of Al ea ; and the function P contains one more free parameter than Pcteqz- 
Pcteq2{x; A 3 , A A ) = 1 + A 3 x Ai ; Pmrs{x; A 3 , A 4 ) = 1 + A 3 yfx + A 4 x. Since two extra de- 
grees of freedom are added, we shall refer to this class of parametrization as (m+2) - i.e. 
minimal plus two. The more general parametrization clearly allows a wider range of varia- 
tion of G(x,Qo). Some pertinent questions are: whether these general parametrizations are 
required by current data; and do these parametrizations give a good indication of the range 
of variation of G(x, Qo)? We shall investigate these questions in some detail in the next two 
sections. 

Finally, although the PDF's determined from global analysis should, in principle, be 
universal, they could, in practice, depend on the choice of data sets - in particular, on the 
choice of u Q cu t" values that specify the minimum hard physical scale (Q, P t , ...) required for 
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data points in the various physical processes to be included in the fit. If the NLO QCD 
theory is truly applicable in the kinematic range of the data, the parton distributions should 
be insensitive to the value of Q cu t- Since current theory does not predict what value Q cut 
should take for each process, this point has to be investigated phenomenologically. 

3 Impact of recent DIS data on the global analysis of 
parton distributions 

Since the publication of the CTEQ3 analysis, more accurate and extensive DIS data from 
NMC @ and HERA |TJ|, as well as new data from E665 have become available. 
In comparing the new data with NLO QCD F 2 computed from CTEQ3M distributions, we 
find general agreement, except for the small-x region where the more precise recent data 
show deviations from the theory curves. This is shown in Fig. [1] for the NMC and HI 
data sets respectively.^] Thus, we first update the CTEQ3 analysis under several different 
conditions, in order to study the impact of these new DIS data on the global analysis of 
parton distributions, especially the extraction of the gluon distribution. 

The magnitude of the uncertainty in G(x, Qq) due to the current uncertainty on a s , will 
be investigated by systematically varying the value of a s over the interval 0.105 < a s < 0.122, 
as mentioned in the previous section. We shall use the short-hand a s for a s (Mz) throughout. 
In terms of QCD Lambda values, this range of a s corresponds to 100 < Af 5 < 280 (MeV) 
and 155 < Af^ < 395 (MeV). We shall in general use the ~MS scheme in NLO QCD. 

To provide a base-line for comparison, we first obtain a series of such fits under identical 
conditions and using the same data sets (i.e. pre-1995) as in the CTEQ3 analysis |13|]. We 
shall refer to this as the A-series.0 By definition, the best fit in this series is the published 
CTEQ3M fit with a s = 0.112 (Alp = 158 MeV). A comparison of the gluon distributions 
that correspond to these values of a s are presented in Fig. || In order to render the differences 
in the various regions of x visible over the range 10~ 4 < x < 1, part (a) highlights the small- 
x region by plotting xG(x,Q) against log x, part (b) accentuates the medium-x range by 
plotting x 2 G(x, Q) vs. log xj 5 ] and part (c) emphasizes large- x by plotting x 2 G(x, Q) vs. x. 
For the many detailed comparisons to follow, these separate plots, though conventional, will 
prove to be rather cumbersome. We consolidate them into one single less-conventional plot 

2 Results are similar for E665 and ZEUS. Comparison with the full data sets will be presented later, cf. 
Fig. | 

3 This series of fits were originally obtained in 1994. They have been used in various phenomenological 
studies related to gluon distributions and a s determination conducted by CTEQ, CDF, and DO Collabora- 
tions. They have not been formally published. 

4 Note that, since x 2 G(x) ■ (flog a; = xG(x) ■ dx = momentum fraction carried within dx, each curve in this 
plot directly depicts the distribution of the momentum fraction carried by the gluon for that set. 



4 



in Fig. ^ in which all curves are normalized by the function x _L5 (l — x) 3 , which takes out 
most of the singular (rapidly vanishing) factors at small (large) x. The scale for the abscissa 
is chosen to be a function of x which smoothly interpolates between log x (at small x) and 
x (at large x) so that the behavior of G(x, Q) over the full x range is more evenly displayed. 
We see that all of the features seen in the three plots of Fig. |3| are evident in this single 
figure. This will be the format of choice in most subsequent comparisons. 

We see in Fig. f| that, in the region x > 0.05 where the largest concentration of data 
used for the fit lie, increasing values of a s lead to decreasing values of G(x, Q) - as expected 
(particularly for the direct photon data) since the product of the two enters into most 
cross-section and evolution kernel formulas.0 As noted before, in the CTEQ3 analysis, and 
therefore in this series of fits, the initial gluon distribution function is parametrized minimally 
as 

G(x, Q ) = A x M (1 - x) M (1 + A 3 x) (2) 

with A\ set to be the same as that of the sea quarks. Hence there are 3 free gluon parameters 
~ ^0,2,3 ~ in the fit. For each a s , we found the best solution to be quite stable against per- 
turbations in the fitting procedure and starting parameters, indicating the parametrization 
and the experimental constraints are well-matched. This also results in an orderly variation 
of G(x, Q) as a s is varied, as seen in the figure. If one takes the range of a s used here as 
representing the current uncertainty on a s , then the spread of the gluon distribution shown 
in Fig. |4] gives the corresponding uncertainty on G(x, Q) (based on the data available prior 
to 1995, and on the variation of a s alone). We should mention that, although quark dis- 
tributions are allowed to vary freely, the valence quark distributions remain practically the 
same for all of the fits in this series, as they are very much pinned down by the precision 
DIS data in the region where they dominate the structure functions. On the other hand, the 
sea quark distributions couple to G(x, Q); thus they do show a systematic variation with a s , 
although the variation is somewhat reduced compared to that of the gluon. 

Next, we investigate the impact of the new DIS data from NMC |18|| , E665 [|ST| and HERA 
ISl on F 2 by repeating the same study, with the new data sets replacing the original ones. 



The resulting series of fits is called the B-series. The quality of these fits (measured in x 2 
values) are similar to those of the A-series. Six representative gluon distributions in this series 
are shown in Fig. |5| along with that of CTEQ3M for reference. It is rather striking to note that 
the spread in G(x, Q) observed above in the small-x (< 0.01) region has been practically 
eliminated. This is precisely the region covered by the HERA experiments. In addition, 
the new gluons are shifted down from those of the A-series in the region 0.05 < x < 0.3 
where all three DIS experiments contribute. At first glance, this may appear surprising 
in view of the conventional wisdom that F 2 data are only sensitive to quarks, not gluons. 



3 The order is reversed for small x, because of the momentum sum rule. 
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However, we must realize that, first, in the small-x region G(x, Q) is quite large — typically 
about 20 times bigger than the quark distributions — thus it has a strong influence, directly 
and indirectly, on all physical quantities through the hard cross-section and the evolution 
equation. Moreover, these fits use the minimal parametrization, including the constraint 
Af = Al ea which strongly couples the behavior of G(x, Q) at small-x to that of sea quarks. 
Thus, the much better determined G(x, Q) just reflects the improved accuracy of new data in 
this region. We note also, the large- x behavior of the new series is somewhat different from 
the A-series, even if there are no new data in that region. This must be due to the indirect 
effect of the required changes below x = 0.1, induced by the restrictive functional form 
Eq. 0, and the constraint imposed by the momentum sum rule. We should point out that 
the absolute value of the gluon distribution in the region above x = 0.5 is very small (about 
10~ 3 compared to its value at x — 0.1); thus the significance of the observed differences 
should not be over-emphasized. 

The minimal parametrization for G(x, Qq) used above was originally chosen in the CTEQ3 
analysis for its economy - all data sets included in these global analyses can be reasonably 
well fitted with this form. This does not prove that the true gluon distribution must fall 
within the range shown above; in particular, the true G(x, Qq) may be more complicated 
than can be represented by this parametrization. (For instance, all global analyses find it 
necessary to use one more parameter to describe the valence quarks.) Only experiments 
probing G(x, Q) in a different way can tell whether our results so far are adequate. Before 
turning to such additional input, we can obtain a different estimate of the uncertainty on 
the gluon distribution that is complementary to the width of the "band" shown in Figs. ^-0. 
We adopt the more general "(m+2)" parametrization of G(x, Q ) already used in CTEQ2: 

G(x, Qq) =Aqx a ^I- x) m (1 + A 3 x M ) (3) 

In addition to introducing the new parameter A 4 compared to Eq. |[ the parameter Ai is 
untied from A s ^ a and treated as free. This results in a new series of fits, called the C-series. 

With two more free parameters than in the B-series, one would expect (i) to fit the 
collective data "better" than before and (ii) to find an increased range of variation of the 
gluon distribution. Indeed, the x 2 f° r the fits decreased slightly (by about 10 (/1000 pts.)) 
compared to the corresponding ones in the B-series. The gluon distributions at Q = 5 GeV 
in this series for 6 values of a s is shown in Fig. First, we see that the range of variation of 
G(x, Q) in this series is much wider as compared to that of series B, although both include 
the same improved DIS data. In particular, in the small-x region the very narrow range in 
series B is very much opened up by the freeing of the A\ parameter for the gluon - since now 
q(x, Q) and a s G(x, Q) can vary independently, the measured F 2 (which depends on both) 
no longer constrains each piece tightly as in the B-series. Secondly, we note that the gluon 
distribution does not vary in a systematic manner as the a s value is varied - in contrast 
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to the well-constrained case in series A and B. Further study has indicated that, unlike in 
the other cases, small changes in the fitting process can lead to different solutions for some 
values of a s . This suggests that the fits are not entirely stable; or, in other words, the system 
becomes somewhat under-constrained with the two extra parameters introduced. 

These observations point to the need for more experimental input in order to better 
measure the gluon distribution. We need new data to determine whether the additional 
degrees of freedom associated with Af and are required for the true gluon or whether 
the restricted form used in series B is already sufficient. If Af and A^ are required, these 
new data could help to stabilize the fits found in the C-series and hence shed light on the 
possible range of G(x, Q) allowed. From the discussion given in the introduction, it is clear 
that inclusive jet production data could be used to help resolve these issues, as we will show 
in the next section. To conclude this section, Table H summarizes the above described three 
series of global fits, as well as those including jet data to be discussed next. 

4 Comparison with New Inclusive Jets Cross-section 

For studying the impact of inclusive jet production cross-section, we use the recent measure- 
ment of da/dEt from the CDF |7j and DO || Collaborations. The preliminary data obtained 
in run IB of the Tevatron by the two experiments are shown in Fig. Although data are 
available for 15 GeV< E t < 450 GeV, we will include in our NLO QCD analysis only data 
above 50 GeV because there are a number of potential theoretical and experimental problems 
that may affect the proper comparison between NLO QCD theory and data for lower E t . 
These include (1) scale uncertainty of NLO QCD calculations, which becomes non- negligible 
at low E t (cf. Fig. |5]a); (2) ambiguities in the definition of the "underlying event" coming 
from the proton-antiproton remnants (cf. Fig. ^b); (3) possible problems in the match be- 
tween theoretical and experimental jet definitions, such as fragmentation products outside 
the jet cone; (4) k t broadening of the initial state partons ||]; and (5) non-perturbative cor- 
rections to the theory, which could be of order 1/E t rather than 1/E^ |22[| . All of these affect 
low E t jets much more than high E t jets, as will be illustrated by two examples, one theoret- 
ical and one experimental. Fig. |8]a shows the scale-dependence of the NLO QCD calculation 
as a function of E t : the theoretical inclusive jet cross-section is shown for several choices 
of the renormalization and factorization scale (// = {Ir — yip) normalized to our standard 
choice \i = E t /2.f\ For low E t , the ratio becomes large and unstable; above 50 — 75 GeV, 
the different choices are within 10% and stay relatively constant — they amount to shifts in 
the overall normalization of the cross-section. Fig. [8|b shows the percentage effect on the 
inclusive jet cross-section in the CDF experiment due to a ±30 % change in the underlying 

6 The theoretical calculations of jet cross-section in this paper are carried out using the EKS program 0. 
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event correction (in Run IA). Again, the uncertainty becomes large below 50 — 75 GeV. 

To emphasize the quantitative aspects of the subsequent analysis, the measured steeply 
falling da/dEt is normalized to the NLO QCD theoretical expectation using the CTEQ3M 
parton distributions (solid horizontal line) and displayed in Fig. || on a linear plot (with 
statistical errors only on the data points). In Fig. ^|, we have taken into account the slightly 
different pseudo-rapidity coverage of the two experiments (0.1 < \rj\ < 0.7 for CDF vs. 
\r)\ < 0.5 for DO) by normalizing each data set with respect to the theory values computed 
with the corresponding 77 range. In addition, we have allowed a small overall normalization of 
the two data sets, well-within the quoted uncertainties, for this comparison. This figure shows 
that the two data sets agree quite well over the entire E t range, especially when considering 
the quoted systematic uncertainties (not shown). See Ref. |23j for more discussions. We will 
discuss the experimental systematic uncertainties in the proper context of the "range" of 
gluon distributions later in this paper. The rise of the data points at high E t values over the 
CTEQ3M expectation, more noticeable for the CDF points, has been the subject of much 
recent discussion and speculation |f7|, |B|, p4fl. We will comment on this issue in the context of 
the global analysis conducted in this paper in a later section. 

Since most inclusive jet data are collected in the central rapidity region, the x-value 
of the PDF's probed is around Xt = 1E t j\fs. For 50 GeV < E t < 450 GeV, the x range 
is approximately 0.06 — 0.5. Over this range, the relative importance of the three parton 
subprocesses - quark-quark, quark-gluon, and gluon-gluon - shifts continuously from being 
gluon-dominated to quark-dominated, as illustrated in Fig. [10]. We should also keep in mind 
that these jet data probe hadron structure at much higher momentum scales than fixed-target 
experiments. Due to the nature of the QCD evolution equation, parton distributions at these 
high momentum scales are determined by those at lower scales and higher x values. Thus the 
effective x-range in G(x, Q ) for some Q , say 1.6 GeV used in our analysis, probed by these 
jet data extends to much higher values than the nominal values mentioned above. Since the 
quark distributions throughout this range are very well pinned down by DIS experiments, 
one expects the jet data to be particularly useful in constraining the gluon distribution. 
The value of a s has considerable influence on the gluon determination for several reasons. 
First, the cross-section for medium x t is proportional to ot 2 s G n (x ) Q) (n = 2,1,0), so that 
as a s increases, G(x, Q) will decrease. Second, a s controls the rate of evolution of G(x, Q) 
and hence affects the slope of the gluon distribution for given measured jet cross-sections. 
Third, a s (/i) itself depends on x through /1 = E t /2 = x^Js/A (at 77 = 0), so that the rate of 
variation of a s (controlled by its strength) is coupled to the x-dependence of G(x, Q) in the 
cross-section formula. 

We now apply the results obtained in Sec. |3| to these jet data to see how the latter agree 
with the predictions of perturbative QCD using these new parton distributions determined 
by the other processes. Fig. [TT] compares the predictions of the PDF's from the B-series 
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(which incorporate the most recent DIS data and use the minimal parameters for the gluon) 
with the jet data, using the same "(Data - Theory) / Theory" format as Fig. |9|. We use the set 
with a s (Mz) = 0.116 as the "Theory" (horizontal solid line) against which the data points 
as well as the predictions of the other fits with different a s values in the series are displayed 
in this plot. To make these comparisons, we allow an overall relative normalization between 
theory and data.[| The normalization factor for the CDF/DO data set ranges from 0.94/0.92 
to 1.08/1.06 for a s = 0.110 to 0.122. |25| The normalization uncertainties quoted by the CDF 
and DO experiments are around 5%. Considering the 7 orders of magnitude of variation of 
the cross-section (Fig. |7]), this is quite remarkable. Within the minimal parametrization of 
the gluon used by the B-series, the parton distributions narrowed down by recent precise DIS 
data (see previous section) are remarkably consistent with the new high statistics inclusive 
hadron-hadron jet data. We also found that the more generally parametrized C-series PDF's 
give qualitatively similar predictions for jet cross-sections compared to the B-series displayed 
in Fig. [11; hence they will not be separately shown. 

The important questions at this point are the following: (i) At a more quantitative level, 
how can these parton distributions be improved by including the jet data in the analysis 
from the beginning; and (ii) will the addition of the jet data reduce the variation of G(x, Q) 
when we use the more general (m+2) parametrization? 



5 New CTEQ parton distribution sets - CTEQ4 

To answer these questions, we have performed an extensive study of the interplay of the 
inclusive jet data with the high-precision DIS and other data within the CTEQ QCD global 
analysis program. The complete set of processes and experiments used is given in Table |^. To 
display explicitly the wide coverage of these experiments over the kinematical variables, we 
show in Fig. |T2]a map of the (x, Q) plane with the data range of the various experiments. We 
see the greatly expanded kinematic coverage compared to a few years ago: in the direction 
of small- x due to the HERA experiments, and in the high Q direction due to the Tevatron 
inclusive jet experiments.^ As before, all processes are treated consistently to NLO accuracy 
in pQCD. This new round of global analysis will be referred to as the CTEQ4 analysis. 

Building upon studies described in the previous sections, we explored all the issues de- 
scribed in Sec. |2|, now with jet data also playing a role. Although the quark distributions 
are coupled to G(x, Q) and a s , they remain tightly constrained by the DIS experiments, 
hence stay very close to those determined before. Thus, our studies concern again mainly 

7 Such a rcnormalization, within errors, is usually allowed in global fitting. 

8 Since these experiments are the only ones in the respective kinematic region, new information on par- 
ton distributions extracted from these data provide challenges to QCD theory for future comparison with 
independent measurements based on other processes. 
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the range of variation of G(x, Q) due to uncertainties in a s and the parametrization of the 
non-perturbative initial distribution. (We have also looked into the influence due to the 
choice of ll Q cu t \ which will be described in the Appendix). Since the results from Sec. [| 
indicate that it is possible to obtain good fits to all the data using the minimal parametriza- 
tion of the gluon distribution even without taking into account the experimental systematic 
errors on the inclusive jet data, we anticipate the most important role of the latter in the 
new analysis is to constrain the possible range of G(x,Q). Hence, we shall use the more 
general (m+2) parametrization which allows a wider range of variation of G(x, Q). We shall 
not include the correlated systematic uncertainties on the jet data since they are not crucial 
for the present purposes. This point will come up again later. More discussions on the 
experimental systematic uncertainties can be found in the Appendix. 

The new generation of CTEQ4 parton distributions are summarized in Table |[ They 
will be described in turn in the following. 



Standard CTEQ4M parton distributions 



We first present the standard fit in the MS scheme which we will designate as the CTEQ4M 
set of parton distributions. The a s {mz) value for this set is 0.116, corresponding to second 
order A5 = 0.202 or A4 = 0.296 GeV. This set gives excellent fit to all data sets. The 
total x 2 f° r 1297 DIS and DY data points is 1320. Detailed information on the x 2 ' s f° r 
the various experiments, in comparison to those obtained using other current and previous 
generations of parton distributions are presented in Tables f| and [| respectively. The direct 
photon and jet data sets are not included in the x 2 table since, without including the sizable 
theoretical uncertainties for the formerf] and experimental systematic errors for the latter, 
the significance of such x 2 values would be difficult to evaluate. The comparison of the CDF 
and DO jet data to the NLO QCD inclusive jet cross-section calculated with the CTEQ4M 
distributions is shown in Fig. [13|. And the comparison of the recent NMC, HI, and ZEUS 
data sets to the fit is shown in Figs. @. 

From Table f|, we see that the CTEQ4M PDF set has the best overall quantitative 
agreement between NLO QCD theory and global data on high energy scattering. It also 
represents a significant improvement over the previous generation of parton distributions, as 
a comparison to Table |5] makes clear. Most of the difference is caused by the new precision 
data from the HERA experiments. Fig. O shows good general agreement of CTEQ4M with 



the jet data, while the much discussed "high E t excess" is still noticeable. We will return to 
this issue in Sec. ^ where an alternative "high E t jet-fit" CTEQ4HJ (included in Table |j) 
will be discussed. Figs. |^ explicitly shows the improvement of CTEQ4M over CTEQ3M in 
describing the recent high-precision DIS experiments. In the Appendix, we will give detailed 



'See the Introduction and Refs. ||, |) for discussions on these uncertainties. 
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information on the parameters which characterize the initial parton distributions at Qo = 1.6 
GeV (which coincides with our choice of the charm threshold). Here, we only note that the 
(Ai,A 2 ) parameters (cf. Eqs. p] & of the gluon and the sea quarks are (—1.21,4.67) and 
(-1.14,8.04) respectively. 



CTEQ4A-series of parton distributions with varying a s and G(x, Q) 

In exploring the range of variation of allowed G(x, Q) by varying the values of a s , changing 
the number of parameters for the gluon, and altering the Q cut of data selection, we have 
found the largest effect is due to the varying of a s . Hence, in presenting a series of PDF's 
which give a reasonable representation of the range of possibilities, we use those generated 
with an a s range centered around the CTEQ4M value of 0.116, which is close to the cur- 
rent world average [|14|. This series will be designated as CTEQ4A-series (shorthand for 
CTEQ4Alpha) — CTEQ4A1 , ...,CTEQ4A5, with CTEQ4A3 being the same as CTEQ4M. 
The x 2 P er point for the 1297 non-jet data points are (1.07,1.02,1.02,1.07,1.19) respectively. 
The higher x 2 values at low values of a s mainly come from the HERA DIS experiments; the 
higher x 2 values at high values of a s are mainly due to the fixed-target DIS experiments |25j . 
The difference in x 2 above minimum, especially for the highest value of a s , is larger than in 
previous CTEQ analyses (e.g. CTEQ2ML vs. CTEQ2M) due to the sharply reduced errors 
on recent DIS data. However, the difference is comparable to that between the MRSJ and 
CTEQ4M x 2 s, cf, Table |j. Because correlations in the experimental errors are not available 
for all experiments, hence have not been included in current global analyses, and since theo- 
retical uncertainties are even harder to quantify, pragmatically, we take these x 2 differences 
as being acceptable for present purposes. 

Fig. [14| shows the comparison of the CTEQ4A parton distribution sets with the two jet 
data sets, using CTEQ4M as the common calibration. The overall normalization factor on 
the jet data sets applied to the various fits range from 0.96 (for CTEQ4A1 on DO points) 
to 1.02 (for CTEQ4A5 on CDF points), well within the experimental uncertainty of ~ 5%. 
Comparing to Fig. [TT] and the range of normalization factors needed there (0.92 — 1.08, which 
is wider than the experimental error), we see the expected improvement of the agreement 
with the jet data. 

The gluon distributions associated with the various values of a s in this series are shown 
in Fig. [T5|. Comparing the CTEQ4A-series to the C-series (same parametrization form for 
G(x, Qo) ), we see that the constraining influence of the jet data has a rather dramatic effect. 
The unstable behavior of the various curves observed in the C-series has been replaced by 
an orderly variation as one steps through the values of a s within the range explored. We 
found, indeed, that for each value of a s , the solution of G(x,Q) is rather unique against 
perturbations in the fitting procedure. 
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One concern is that the variation in the CTEQ4A series is too small due to the lack of 
treatment of systematic uncertainties in the jet data. To address this issue, we compare 
the change in the calculated jet cross-sections between the extremes of the CTEQ4A series 



to the largest .E^-dependent uncertainty in the CDF data. See Fig. [16]. It shows that the 
range of variation in the CTEQ4A series is about 10% in the moderate Et range, while 
the experimental systematic uncertainty is about the same-P^ This observation lends some 
confidence that this series gives a reasonable estimate of the range of variation of G(x,Q). 
To the extent that there are sources of uncertainty other than a s , the variation in G(x, Q) 
given here may be considered a minimum range. However, our study does indicate that the 
variation due to the uncertainty of a s may be the dominant one. 

Fig. [17] shows a comparison of some of the new gluon and singlet quark distributions with 
those of CTEQ3M and MRSJ in the usual form xf(x,Q) without the normalization factor 
as in previous figures. On this conventional plot, differences in G(x, Q) can be seen only in 
the small- x region, and the CTEQ3M and CTEQ4M gluons appear to be indistinguishable. 
Differences in the singlet quark distribution are more evident near x = 0.01. The fact that 
only small changes in the parton distributions result from adding so much new data in the 
global analysis is testament to the impressive progress in pinning down these parton distri- 
butions that has been made in recent years. These changes, though small, are nonetheless 
physically significant, as demonstrated by the substantial differences in y 2 values between 
the new and old parton distribution sets on the precision experiments given in Tables f| 
and ||. 



Other CTEQ4 parton distributions 

Along with the standard CTEQ4M MS parton distributions, we have also obtained corre- 
sponding parton distributions in the "DIS scheme" — CTEQ4D. CTEQ4D uses the same 
value of a s (= 0.116) as CTEQ4M; it is obtained by fitting under identical conditions as 
CTEQ4M except that the hard cross-sections are evaluated in the DIS scheme. The \ 2 
values of this fit are comparable to those of CTEQ4M. In addition to these two standard 
sets, for applications requiring leading order (LO) calculations and low values of the scale Q 
(LQ), we also provide appropriate parton distribution sets labelled CTEQ4L and CTEQ4LQ 
respectively. The CTEQ4LQ set can be used for Q 2 > Q\ = 0.5 GeV 2 . Q It was obtained by 
fitting the same data sets as the other PDF sets. Since the proper treatment of low Q data 
must involve more physics input (such as higher twist effects) than included here, CTEQ4LQ 

10 Of course, given the good agreement between the two Tevatron experiments p3| ], if the CDF jet data 
requires a significant change due to a systematic error, the DO data would require the same change, an 
unlikely occurrence since there is almost no correlation in the two experimental measurements. 

11 Bclow this scale (Q = 0.7 GeV) the QCD coupling approaches unity, the perturbative formulas certainly 
cease to be meaningful. 
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represents only an extrapolation of twist-two QCD physics into the low Q region — it is not 
intended to be a best fit. However, as demonstrated by the GRV parton distribution sets 
|P2|| , this kind of extrapolation often turns out to compare rather well with data in the low 
Q region. Comparison of CTEQ4LQ structure functions to the NMC, E665 and HI data in 
the range 1.0 < Q < 3.0 GeV is shown in Fig. ||. The parameters for CTEQ4D, CTEQ4L 
and CTEQ4LQ are also given in the Appendix. 

The rather remarkably consistent picture resulting from this round of CTEQ4 global 
analysis incorporating jet data from hadron collisions provides a new generation of improved 
parton distributions for making calculations and predictions on high energy processes both 
within and beyond the standard model. The more tightly constrained parton distributions 
can also lay the foundation for more stringent tests of the pQCD framework and provide the 
basis for discerning signals of new physics. 

At present, a remaining area of some uncertainty is the gluon distribution in the "large 
x" region, beyond say 0.25, where neither the DIS nor the direct photon data give tight 
constraints. For the DIS process, the sensitivity to the gluon begins below x = 0.1. For the 
direct photon process there are a number of theoretical uncertainties which are not yet under 
control, as already discussed in the Introduction. The noticeable rise of the inclusive jet data 
points J7[ above all "theory" curves shown so far may be related to the conventional choices 
of parametrization of the non-perturbative function G(x,Qi), which restricts its behavior 
in the large x region. This possibility, first raised in Ref. M, will be discussed next in the 
context of the CTEQ4 analysis presented above. 



6 High E t Jets and Parton Distributions 

The higher-than-expected inclusive jet cross-sections, first measured by the CDF collabora- 
tion [7|] for E t > 200 GeV, were observed in comparison to the existing parton distribution 
sets, including CTEQ3M as shown in Fig. [| This "excess" is reduced slightly when jet data 
are included in the global fit, but is still noticeable in Figs. [TBI and IT3 for the CTEQ4A series 



of distributions. This is understandable since the high E t data points have large errors, so do 
not carry much statistical weight in the fitting process, and the simple (unsigned) x 2 is n °t 
sensitive to the observed pattern that all the points are higher than the theoretical prediction 
in the large E t region. Ref. || investigated the feasibility of accommodating these higher 
cross-sections in the conventional QCD framework by exploiting the flexibility of G(x, Q) 
at higher values of x where there are few independent constraints, while maintaining the 
agreement with other data sets in the global analysis. To do this, it is necessary to (i) 
provide enough flexibility in the parametrization of G(x, Qq) to allow for behaviors different 
from the usual (but arbitrary) choice; and (ii) focus on the high E t data points and assign 
them more statistical weight than their nominal values in order to force a better agreement 
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between theory and experiment. Thus, the spirit of the investigation is not to obtain a "best 
fit" in the usual sense. Rather, it is (i) to find out whether such solutions exist; and (ii) if 
they do exist, to quantify how well these solutions agree with other data sets as compared to 
conventional parton distribution sets. The global analysis work described in Sec. |5] without 
special attention to the high E t points provides the natural setting to put the results of 
Ref. H in context. 

Ref. was performed using the CDF Run-IA data — the only high statistics inclusive 
jet measurement available at the time. Two illustrative "solutions" of the type described 
above were presented — one with the normalization fixed at 1.0 with respect to the CDF 
data, the other with a normalization factor of 0.93. Fig. [19] compares predictions of the 
normalization=1.0 PDF set, which we shall refer to as the CTEQ4HJ set, with the more 
recent Run-IB results of both CDF and DO. For this comparison, an overall normalization 
factor of 1.01(0.98) for the CDF(DO) data set is found to be optimal in bringing agreement 
between theory and experiment .0 The consistency between the two data sets, as well as be- 
tween theory and experiment, displayed by this comparison appears to be rather remarkable 
(again, bearing in mind the neglect of systematic errors other than overall normalization). 
Results shown in Table |4] quantify the \ 2 values obtained while accommodating the high 
E t jets in the global fit in this particular case. Compared to the best fit CTEQ4M, the 
overall \ 2 f° r CTEQ4HJ is indeed slightly higher. But this difference is much smaller than 
the differences discussed earlier in the CTEQ4A series, and much smaller than the difference 
between MRSJ and CTEQ4M. Thus the price for accommodating the high E t jets is negligi- 
ble. In addition, the difference between CTEQ4HJ and CTEQ4M is almost entirely due to 
the BCDMS data, even though the BCDMS x 2 for CTEQ4HJ by itself is quite good. This 
change is due to the fact that, in the CTEQ4M fit, the BCDMS data set is the dominant 
one determining the large- x quark distributions, while, in the CTEQ4HJ fit, the jet data 
set is in competition for these quark parameters, and they are changed by minute amounts. 
This is shown in Fig. ^ where the residuals between BCDMS data and theory are shown for 
CTEQ4M and CTEQ4HJ. The residuals are almost identical, which, together with Table [|, 
confirms the fact that even though CTEQ4HJ does not give the absolute overall best fit to 
all data, it provides an extremely good description of all data sets. It should be considered 
as a candidate for the gluon distribution in nature.^ In the future we will need strong, 
independent measurements of the large-x gluons in order to clarify the situation with the 
high-i?t jets. 

12 The change of CDF normalization factor from 1.0 to 1.01 is attributable to the switch from the Run-IA 
to the Run-IB data set. 

13 This is to be contrasted with the conclusion of incompatibility between the inclusive jet and DIS data 
reached by Ref. j53). Their fit to inclusive jet data over the full E t range (the MRSJ' set) gives rise to an 
extremely large \ 2 f° r the BCDMS data set. 
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7 Summary 



In this study of the impact of recent DIS and inclusive jet data on the global QCD analysis of 
lepton-hadron and hadron-hadron processes, we see significant progress in demonstrating the 
consistency of the NLO QCD framework, and in narrowing the uncertainties on the elusive 
but important gluon distribution. Specifically, 

• The recent NMC, E665, HI and ZEUS data considerably narrow down parton distri- 
butions and limit the behavior of the gluon, especially if one uses the minimal form of 
the gluon parameterization used by CTEQ3; 

• The new inclusive jet data agree well with theory predictions based on PDF's deter- 
mined by the other processes, with the possible exception of the high E t data points. 

• By adding jet data to the global analysis, it is possible to further explore the range of 
variation of the gluon distribution using a more general parametrization. Although the 
jet data set covers a limited x-region, its effect is felt over the entire x-range - because 
it complements the other data sets well. 

• Based on these investigations, a new generation of CTEQ4 parton distributions for a 
variety of features are presented: they are tabulated in Table ^. 

• Three sources contributing to the uncertainty of the gluon distribution have been in- 
vestigated: (i) by letting a s vary over its current range of uncertainty; (ii) by increasing 
the degree of freedom for parametrizing the non-perturbative initial gluon distribution, 
and (iii) by varying the Q cu t in selecting data for the global fits. The largest effect is 
due to a s . 

• These studies help to delineate the range of variation of G(x, Q) over the range 10~ 4 < 
x < 0.25. Further work is needed in exploring the range of uncertainty of the gluon 
and other parton distributions by systematically varying the relevant parameters of 
the global analysis. 

• For larger values of x, more definitive experimental results on inclusive jet and direct 
photon production as well as improved theory are needed for further progress. The 
observed high p t "excess" jet cross-section can be accommodated by a modified gluon 
distribution, represented by the CTEQ4HJ set, since no other independent measure- 
ment constrains it in this range. 

In view of the strong correlation between the gluon distribution and a s , narrowing the 
uncertainty in the latter will significantly improve the determination of G(x, Q). What can a 
global analysis of experimental data described in this paper contribute to the measurement 
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of a s l To explore this question, one needs to study in some detail the sensitivity of each 
process which contributes to the global analysis to the variation of a s . This problem will be 
pursued in a separate analysis. 



Appendix 

CTEQ4 Parton Distribution Parameters 

The initial parton distributions at Q — Qq, f l (x, Qq), are parametrized in general as in Eq. ^| 
for the gluon G and the quark flavors d v , u v , u + d, s (s); except for the combination d — u 
(which does not have to be positive definite) which is parametrized as: 



d-u = A x Al (1 -x) M (1 + A 3 y/x + Ai 



x 



For all parton distribution sets, Qq = 1.6 GeV, except for CTEQ4LQ which has Qo = 0.7 
GeV. Tables of the coefficients {A l n ; n = 1, ..,4; i = flavors} for the three standard parton 
distribution sets CTEQ4M, CTEQ4D, CTEQ4L and the low-Q set CTEQ4LQ are given 
below, in Tables |6|,|7], ^| and All parton distribution sets listed in Table [| are available in 
fortran program form by request F] or via WWW at [http:/ /www.phys.psu.edu/~cteq/. 



Experimental Normalization Factors 

The x 2 tables HJ5| are obtained by allowing the experimental data sets to "float" with re- 
spect to the theory cross-sections. For CTEQ distributions, a x 2 penalty is included in 
the fitting process for deviations of the normalization factors with respect to the respective 
overall experimental normalization errors. For non-CTEQ distributions, we simply obtained 
the minimum x 2 by varying the normalization factors without such penalty. The resulting 
normalization factors which go with Tables are given in Tables [TO] and |TT| . 

Experimental Systematic Uncertainties 

For DIS, DY and direct photon data, we follow the usual procedure of combining in quadra- 
ture point-to-point systematic errors given by the experiments with the statistical errors. 
Correlated systematic errors other than overall normalization are not generally available 
from most experiments. For a few where they are, we have done separate studies of the 
consequences of incorporating them in the global analysis and found they do not affect the 
best fit parameters by any significant amount. See Ref. fl3f . 

For the preliminary inclusive jet data, only the normalization uncertainty is taken into 
account in the global fit. The rationale has been explained in Sec. |5|. The fully correlated 

14 Requests can be sent to Lai_H@Pa.Msu.Edu or Tung@Pa.Msu.Edu. 
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systematic errors from CDF, although available, are not easily implemented in a way which is 
consistent with all the other data sets. (A separate study on the effects of these uncertainties 
employing the full correlation matrix is underway, and it will be reported in the future.) 
These errors are not yet available for the DO data set. The main effect of omitting the 
systematic errors on jets is to increase somewhat the relative weight of this data set in the 
global analysis. This will not affect the fits substantially because the jet data agree well with 
parton distributions determined from other processes, as discussed in Sec. f|. 

In general, the question of assigning appropriate relative weights to different experimen- 
tal data sets in a global analysis is a difficult one. An experiment with few data points 
which is however particularly sensitive to some physical parameters than all the others can 
sometimes be emphasized justifiably in a global x 2 minimization process, otherwise it will 
be overwhelmed by the far more numerous data sets and the sensitivity will be lost. As 
an extreme example, the NA51 experiment [p9fl , which has an important impact on the de- 
termination of the flavor SU(2) assymmetry of the sea quarks (u — d), consists of only one 
data point. It has to be appropriately emphasized in a global analysis to have an effect in 
differentiating the sea quarks. 



Dependence on the Choice of Q 



cut 



In all global QCD studies, a set of cut-offs on the hard scale "Q" for various processes is used 
in data selection. In recent CTEQ analyses, this Q cu t has been 2 GeV on Q and 3.5 GeV 
on W in DIS, 2 GeV on Q (the invariant lepton pair mass) in Drell-Yan process, and 4 GeV 
on p t in direct photon production. As a final check on the reliability of the results described 
in the previous section, we test the sensitivity of the fits to the value of these cut-offs in 
order to gauge possible influence due to non-perturbative or higher-twist effects. For this 
purpose, we carried out several series of analyses similar to the CTEQ4A-series above, but 
with the minimum Q cut raised progressively from 2 GeV to 3, 4, and 5 GeV. Data points 
excluded by these higher Q cut's are mainly those of fixed-target DIS experiments. We found 



our results to be rather stable under these changes. Fig. |2l] compares the gluon distributions 
from three PDF sets obtained with three Q cut values mentioned above, all for a given a s 
value of 0.113. We see that the differences are quite small - smaller than those due to the 
variation of a s (with the same Q cu t) shown in Fig. [TJ] and described in Sec. |5[ The subtle 



differences, especially in relation to sensitivity on a s values, will be discussed elsewhere JTB| . 

15 This issue has previously been investigated in Ref. |?4j. The accuracy of both experiments and theory 
have improved dramatically since then. 
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Figure 1 : Comparison of NLO calculations based on the previous generation CTEQ3M 
parton distributions with the latest NMC (a) and HI (b) data in the small-x region where 
discrepencies appear. 



21 




g 2 4 6 8 10 1 10 1 

& Q(GeV) Q(GeV) Q(GeV) 

o' 


cn 




10" 4 10" 3 10" 2 10" 1 



X 

0.15 I 




0.2 0.4 0.6 0.8 

X 



Figure 3 : Series-A gluon distributions in the small-, medium-, and large- £ regions. A. 105 
refers to the gluon associated with a s (M z ) = 0.105, and likewise for the other ones. 
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to display clearly the behavior of G(x, Q) over the entire x-range. For the same purpose, the 
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Figure 5 : Series-B gluon distributions normalized by the function x L5 (l — x) 3 (cf. caption 
of previous figure.) 
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Figure 6 : Series-C gluon distributions normalized by the function x L5 (l — x) c 
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Figure 7 : Inclusive jet cross-section measured by the CDF and DO collaborations in Run- 
IB at the Tevatron. (Averaged over 0.1 < \rj\ < 0.7 in the case of CDF and \rj\ < 0.5 in the 
case of DO.) 
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Figure 8 : Two examples of sources of uncertainties in comparing inclusive jet data with 
NLO QCD theory: (a) Fractional difference between da(E t , fi) / dE t and da(E t , fi = E t /2)/ dE t 
(for the CDF rapidity coverage 0.1 < < 0.7) as a function of E t for a variety values of 
(b) Fractional change in the cross-section due to ±30 % change in underlying event correction 
in the CDF experiment. 
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Figure 9 : Inclusive jet cross-section measured by the CDF and DO collaborations in Run- 
IB at the Tevatron normalized to NLO QCD calculations based on CTEQ3M PDF's. The 
difference in rapidity coverage of the two experiments is taken into account. 
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Figure 10 : Relative contribution to the inclusive jet cross-section due to the various par- 
tonic subprocesses. 
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Figure 11 : Inclusive jet cross-section of CDF and DO compared to NLO QCD calculations 
based on the new B-series parton distributions. 
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Figure 12 : Kinematic map in the (x, Q) plane of data points used in the current global 
analysis. 
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Figure 13 : Inclusive jet cross-section of CDF and DO compared to NLO QCD calculations 
based on the new CTEQ4M parton distributions. 
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Figure 14 : Inclusive jet cross-section of CDF and DO compared to NLO QCD calculations 
based on the new CTEQ4A series of parton distributions. 




Figure 15 : Series-CTEQ4A gluon distributions normalized by the function x L5 (l — x) 3 . 
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Figure 16 : Percentage range of variation of the inclusive jet cross-section from the two 
extreme CTEQ4A PDF sets (CTEQ4A1 and CTEQ4A5) compared to the largest of the E t 
dependent systematic uncertainties. 




Figure 17 : Comparison of xG(x,Q) and x S(x,Q) between some new parton distribution 
sets and those from CTEQ3M. S(x, Q) is the singlet quark distribution (sum over all flavors). 
The CTEQ3M and CTEQ4M gluons appear to lie on top of each other. The same is true 
for the CTEQ4A1 and MRSJ gluons. Differences in G(x, Q) for x > 0.01 are not evident in 
this plot. 
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Figure 19 : Inclusive jet cross-section of CDF and DO compared to NLO QCD calculations 
based on the CTEQ4HJ parton distributions. 
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Figure 20 : Percentage deviation of BCDMS proton data from NLO QCD values based on 
CTEQ4M and CTEQ4HJ. Both PDF sets give good fits. 
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Figure 21 : Comparison of gluon distributions obtained in three global fits using three 
different values of Q cut in data selection. 
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Section 




DIS data 


Jet Data 


trization 


discussed 


A 
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m 


Si 


C 


X 




m+2 




CTEQ4A 
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Qcut 


X 


X 
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Appendix 



Table 1: Several series of global fits on which the physics discussions are based. "New DIS 
data" refers to those becoming available since 1995. Minimal parametrization "m" refers to 
Eq. 2; and "m+2" refers to Eq. 3. The last column refers to the section number where the 
specific series is discussed. 
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Process 


Experiment 


Measurable 


Data Points 


Ref. 


DIS 


BCDMS 


r 2 Hi r 2 D 


324 


\m 




NMC 


Z?M 

r 2 Hi r 2 Di r 2 n/p 


297 


|18| 




E665 


r 2 Hi r 2 D 


70 


|21 






HI 


F 2 H 


172 


|1S| 




ZEUS 


F 2 H 


179 


m 




CCFR 


r 2 Fei x r 3 Fe 


126 


\21\ 


Drell-Yan 


E605 


sda/d^/rdy 


119 


|2E| 




NA-51 


A DY 


1 


1 29 




W-prod. 


CDF 


Lepton asym. 


9 


[0 


Direct 7 


WA70 


Ed 3 a/d 3 p 


8 




1 






UA6 


Ed 3 a/d 3 p 


16 






Incl. Jet 


CDF 


da/dE t 


36 










DO 


da/dEt 


26 




§ 



Table 2: List of processes and experiments used in the Global analysis. 



PDF set 


Description 


a s (m z ) 


Ql (GeV 2 ) 




Standard Sets 






CTEQ4M 


MS scheme 


0.116 


2.56 


CTEQ4D 


DIS scheme 


0.116 


2.56 


CTEQ4L 


Leading Order 


0.132 


2.56 




a s series 






CTEQ4A1 


1 


0.110 


2.56 


CTEQ4A2 


2 


0.113 


2.56 


CTEQ4A3 


Same as CTEQ4M 


0.116 


2.56 


CTEQ4A4 


4 


0.119 


2.56 


CTEQ4A5 


5 


0.122 


2.56 




Specials 






CTEQ4HJ 


"Hi- Jet" 


0.116 


2.56 


CTEQ4LQ 


"Low Q " 


0.114 


0.49 



Table 3: List of new CTEQ4 parton distributions and their characteristics. 
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Expt. 


#Dts 


CTEQ4M 


CTEQ4HJ 


CTEQ4LQ 


MRS J 


DM Q H 


lOo 


1/1/1 sfn 

14t4t.ol U.ou ) 


17*? f\(1 CfV\ 
1 I O.U^l.UO ) 




lOQ 1/1 f\Q\ 

loo. i ^i.uy ) 


1 J \_j 1 y iv 1 kJ 


1 5fi 


1 85 fid 1 Q"l 


905 9f1 39"l 




999 3(1 4.7"! 


NMC H 


104 


97.3(0.94) 


91.7(0.88) 


96.0(0.92) 


113.4(1.09) 


NMC D 


104 


93.3(0.90) 


90.2(0.87) 


97.9(0.94) 


122.7(1.18) 


NMC R 


89 


130.8(1.47) 


133.5(1.50) 


132.6(1.49) 


142.4(1.60) 


EQQ5 H 


35 


41.3(1.18) 


38.5(1.10) 


44.5(1.27) 


37.8(1.08) 


E665 D 


35 


32.3(0.92) 


33.5(0.96) 


34.3(0.98) 


29.8(0.85) 


CCFRF 2 


63 


83.2(1.32) 


72.4(1.15) 


74.3(1.18) 


107.7(1.71) 


CCFRF 3 


63 


46.5(0.74) 


45.5(0.72) 


49.9(0.79) 


57.8(0.92) 


ZEUS 


179 


243.4(1.36) 


232.7(1.30) 


268.5(1.50) 


252.4(1.41) 


HI 


172 


118.9(0.69) 


120.2(0.70) 


131.9(0.77) 


109.6(0.64) 


CDFA W 


9 


4.3(0.48) 


3.4(0.38) 


3.8(0.42) 


3.3(0.37) 


NA51 


1 


0.6(0.63) 


0.5(0.49) 


0.4(0.41) 


2.5(2.47) 


£605 


119 


97.7(0.82) 


101.6(0.85) 


100.4(0.84) 


97.8(0.82) 


Total 


1297 


1320 


1343 


1356 


1490 



Table 4: Total x 2 values and their distribution among the DIS and DY experiments for 
current generation of parton distributions which take into account the most recent HERA 
(1996) and NMC (1995) data. In parentheses are the xV^omt values. 
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Expt. 


#pts 


MRSA' 


CTEQ3M 


MRSA 


GRV 


BCDMS H 


168 


156.9(0.93) 


128.7(0.77) 


168.0(1.00) 


250.3(1.49) 


BCDMS D 


156 


213.7(1.37) 


190.3(1.22) 


215.3(1.38) 


187.2(1.20) 


NMC H 


104 


129.0(1.24) 


146.6(1.41) 


114.4(1.10) 


123.8(1.19) 


NMC D 


104 


151.8(1.46) 


137.3(1.32) 


135.2(1.30) 


115.4(1.11) 


NMC/j 


89 


143.3(1.61) 


134.4(1.51) 


140.6(1.58) 


129.0(1.45) 


E665 H 


35 


38.2(1.09) 


47.6(1.36) 


37.8(1.08) 


39.9(1.14) 


E665 D 


35 


29.1(0.83) 


44.5(1.27) 


29.5(0.84) 


29.8(0.85) 


CCFR F 2 


63 


68.0(1.08) 


66.2(1.05) 


68.7(1.09) 


164.4(2.61) 


CCFR F 3 


63 


54.1(0.86) 


41.9(0.67) 


61.7(0.98) 


114.7(1.82) 


ZEUS 


179 


368.7(2.06) 


549.5(3.07) 


1222.6(6.83) 


843.1(4.71) 


HI 


172 


149.5(0.87) 


220.2(1.28) 


407.6(2.37) 


404.2(2.35) 


CDF A w 


9 


4.2(0.47) 


3.0(0.33) 


3.7(0.41) 


9.6(1.07) 


NA51 


1 


0.1(0.06) 


0.4(0.42) 


0.01(0.01) 


0.01(0.01) 


E605 


119 


93.5(0.79) 


92.6(0.78) 


95.9(0.81) 


90.3(0.76) 


Total 


1297 


1600 


1803 


2701 


2502 



Table 5: Total x 2 values and their distribution among the DIS and DY experiments for the 
previous generation of parton distributions which includes experimental data available in 
1995 (MRSA') or before 1995 (CTEQ3M, MRSA). GRV does not perform a full global fit. 
Since it is used widely, it is included here for reference. In parantheses are the \ 2 1 'point 
values. 



Parton 


A 


A l 


A 2 


^3 


A, 


% Momentum 




0.640 


0.501 


4.247 


2.690 


0.333 


11.2 




1.344 


0.501 


3.689 


6.402 


0.873 


30.6 


xg 


1.123 


-0.206 


4.673 


4.269 


1.508 


41.7 


x(d — u) 


0.071 


0.501 


8.041 


0.000 


30.000 




x(d + u) 


0.255 


-0.143 


8.041 


6.112 


1.000 


13.2 


xs 


0.064 


-0.143 


8.041 


6.112 


1.000 


3.3 



Table 6: Parameters for the CTEQ4M initial parton distributions at Q = 1.6 GeV. Also, 
tts(^z) = 0.116, corresponding to A 5 = 202 MeV. 
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Parton 


A 


Ai 


A 2 


^3 


A, 


% Momentum 


xd v 


0.724 


0.490 


3.839 


1.688 


0.338 


11.3 




1.528 


0.490 


3.554 


6.448 


1.162 


30.4 


xg 


2.141 


-0.058 


7.554 


36.405 


2.223 


43.7 


x(d — u) 


0.054 


0.490 


7.200 


0.000 


30.000 




x(d + u) 


0.154 


-0.227 


7.200 


6.949 


1.000 


11.7 


xs 


0.038 


-0.227 


7.200 


6.949 


1.000 


2.9 



Table 7: Parameters for the CTEQ4D initial parton distributions at Q — 1.6 GeV. Also, 
a s (m z ) = 0.116, corresponding to NLO A 5 = 202 MeV. 



Parton 


A 


A 1 


A 2 


^3 


A, 


% Momentum 




0.702 


0.443 


4.003 


2.433 


0.622 


10.9 




1.226 


0.443 


3.465 


7.589 


1.146 


30.1 


xg 


0.854 


-0.305 


3.666 


1.846 


1.968 


41.8 


x(d — u) 


0.050 


0.443 


6.877 


0.000 


30.000 




x(d + u) 


0.201 


-0.200 


6.877 


5.644 


1.000 


13.8 


xs 


0.050 


-0.200 


6.877 


5.644 


1.000 


3.5 



Table 8: Parameters for the CTEQ4L initial parton distributions at Q = 1.6 GeV. Also, 
LO A 5 = 181 MeV 



Parton 


A 


A l 


A 2 


^3 


A A 


% Momentum 


tic d><^j 


0.852 


0.573 


4.060 


4.852 


0.693 


14.7 




1.315 


0.573 


3.281 


10.614 


1.034 


40.4 


xg 


39.873 


1.889 


5.389 


0.618 


0.474 


31.2 


x(d — u) 


0.093 


0.573 


7.293 


0.000 


30.000 




x(d + u) 


0.578 


0.143 


7.293 


1.858 


1.000 


11.7 


xs 


0.096 


0.143 


7.293 


1.858 


1.000 


1.9 



Table 9: Parameters for the CTEQ4LQ initial parton distributions at Q = 0.7 GeV. Also, 
NLO A 5 = 174 MeV 
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Expt. 


CTEQ4M 


CTEQ4HJ 


CTEQ4LQ 


MRSJ 


BCDMS 


0.988 


0.983 


0.993 


0.978 


NMC 


1.016 


1.015 


1.022 


1.018 


NMC R 


1.000 


1.000 


1.000 


1.000 


£665 


1.013 


1.027 


1.009 


1.041 


CCFR 


0.976 


0.971 


0.983 


0.968 


ZEUS 


1.004 


0.999 


1.001 


1.018 


HI 


0.993 


0.978 


0.987 


0.994 


CDFA W 


1.000 


1.000 


1.000 


1.000 


NASI 


1.000 


1.000 


1.000 


1.000 


£605 


1.076 


1.051 


1.075 


1.070 



Table 10: List of normalization factors for the experiments which minimize the x 2 's given in 
the corresponding x 2 table. 



Expt. 


MRSA' 


CTEQ3M 


MRSA 


GRV 


BCDMS 


0.977 


0.988 


0.977 


0.957 


NMC 


1.019 


1.005 


1.018 


0.988 


NMC R 


1.000 


1.000 


1.000 


1.000 


£665 


1.045 


0.997 


1.040 


0.992 


CCFR 


0.968 


0.976 


0.968 


0.949 


ZEUS 


1.023 


0.995 


1.099 


0.878 


HI 


0.988 


0.957 


1.030 


0.836 


CDFA W 


1.000 


1.000 


1.000 


1.000 


NA51 


1.000 


1.000 


1.000 


1.000 


£605 


1.025 


1.097 


1.008 


1.012 



Table 11: List of normalization factors for the experiments which minimize the x 2 's given in 
the corresponding x 2 table. 
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